Pulsed squeezed light: simultaneous squeezing of múltiple modes 
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We analyze the spectral properties of squeezed light produced by means of pulsed, single-pass 
degenerate parametric down-conversion. The multimode output of this process can be decomposed 
into characteristic modes undergoing independent squeezing evolution akin to the Schmidt decompo- 
sition of the biphoton spectrum. The main features of this decomposition can be understood using 
a simple analytical model developed in the perturbative regime. In the strong pumping regime, 
for which the perturbative approach is not vàlid, we present a numerical analysis, specializing to 
the case of one-dimensional propagation in a beta-barium borate waveguide. Characterization of 
the squeezing modes provides us with an insight necessary for optimizing homodyne detection of 
squeezing. For a weak parametric process, efficient squeezing is found in a broad range of local 
oscillator modes, whereas the intense generation regime places much more stringent conditions on 
the local oscillator. We point out that without meeting these conditions, the detected squeezing can 
actually diminish with the increasing pumping strength, and we expose physical reasons behind this 
inefficiency. 

PACS numbers: 42.50.Dv, 03.67.-a 



I. INTRODUCTION 

Squeezed states of the radiation field are those in which 
the field fluctuations in one of the quadratures are re- 
duced below the vacuum noise level. They were among 
the first nonclassical states of light generated experimen- 
tally 0,0 hi the 1980s, and were considered promising as 
a way of overcoming the shot-noise precision restrictions 
in optical measurements 0] and enhancing the capacity 
of communication channels Nowadays, the interest 
to generate squeezed light is mainly due to its applica- 
bility in quantum information and communication. In 
this context it is viewed primarily as a means of gener- 
ating quadrature entanglement a valuable resource 
that can be applied, for example, for quantum telepor- 
tation jfj, computation Q, cryptography 0] and manip- 
ulation of atòmic quantum states There has been 
also a theoretical effort towards developing protocols for 
the continuous-variable error correction 0, EI and en- 
tanglement purification (Ï3 |. Furthermore, operations 
that enable entanglement distillation in the continuous- 
variable domain have been described (ï^ . 

The present paper concentrates on one experimental 
method for producing optical squeezing, namely an opti- 
cal parametric amplifier (OPA), implementing degener- 
ate parametric down-conversion in a single-pass, pulsed 
configuration. A strong pump pulse propagating through 
a nonlinear crystal generates pairs of down-converted 
photons. The idealized experimental arrangement would 
be such that the photons in each generated pair are 
completely indistinguishable and contribute to the same 
spatio-temporal mode of the generated field. A combina- 
tion of a high pump pulse energy and high optical non- 
linearity leads to múltiple pair creation events, which can 
be described in the Heisenberg picture as a linear trans- 



formation between the field operators at the input ài n 
and the output à ou t of the amplifier: 

à ut = à in coshC + à\ n sinhC- (1) 

When the initial state of the field is vacuum, the above 
equation defines a squeezed state with the degree of 
squeezing £, corresponding to the mean square quadra- 
ture noise scaling factor equal to e -2 *». A crucial as- 
sumption underlying Eq. frequently used in theoret- 
ical works, is that all the down-converted photons are 
emitted to a single mode described by the annihilation 
operator à out . However, this does not correspond to a 
realistic experimental situation. 

In an experiment, pulsed squeezed light was realized 
for the first time by Slusher et al. 0] and has since been 
reproduced and improved by a number of researchers 
[ïfll fiti IrR íïïi. lïflj . In all these experiments, a pulsed 
làser is first frequency doubled and then undergoes para- 
metric down-conversion in a spatially- and spectrally- 
degenerate type I configuration. The generated optical 
state is then characterized by means of homodyne de- 
tection. In contrast to the continuous-wave regime, the 
pulsed squeezed mode does not associate with a sideband 
of the carrier wave, but is localized within the tempo- 
ral frame of the pump làser pulse (2ÏJ]. This appealing 
feature is however compromised by a complex spectral 
structure of the generated photon pairs. In single-pass 
down-conversion, biphotons are emitted not necessarily 
into a particular mode, but into a wide range of modes 
restricted only by the energy-conservation and phase- 
matching conditions. The usual absence of a dominant 
spatio-temporal mode can lead to difficulties when we at- 
tempt to observe quantum interference of the generated 
squeezed vacuum with other classical or nonclassical op- 
tical states. The negative effects range from efficiency 
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reduction in homodyne detection to extra dark counts in 
single photon measurements. 

The present paper is dedicated to studying the spectral 
mode of the pulsed squeezed state. A particular question 
we aim to answer is the following: can we deflne a partic- 
ular mode in which high-efficiency squeezing is achieved? 
Are there methods to engineer this mode and what are 
the most efficient ways to access it in homodyne mea- 
surements? The arguments given in the previous para- 
graph may lead one to believe that there is no mode in 
which high-efficiency squeezing is present. Surprisingly, 
we demonstrate just the opposite: we find modes in which 
pure squeezed states are generated, and also a whole con- 
tinuum of their phase-locked superpositions turns out to 
become efficiently squeezed. In fact, good squeezing is 
present in a broad, continuous range of pulsed modes 
that are temporally synchronized with the pump. 

The bàsic theoretical tool in our analysis will be 
the Bloch-Messiah reduction of an arbitrary Bogoliubov 
transformation of bosonic fleld operators [21]. This pro- 
cedure will allow us to identify a set of independent 
squeezing modes whose evolution takes a simple form 
analogous to Eq. QJ. Such a decomposition results di- 
rectly from the preservation of the commutation relations 
between the field operators in course of the parametric 
evolution. This fact allows us to apply the numerical 
technique of singular- value decomposition to integral ker- 
nels describing the transformation of the field operators, 
which provides effective means to analyze a realistic para- 
metric amplifier. The decomposition into independent 
squeezers has been given previously by Bennik and Boyd 
|22l |. who deduced it using intricate operator algebraic 
methods 0. 

As a physical example, we analyse one-dimensional 
propagation in a beta-barium borate (BBO) crystal serv- 
ing as a model of a waveguide amplifier. The nonlin- 
earity of a BBO crystal has an instantaneous character 
[2Í|, as this medium is transparent from infrared to ul- 
traviolet and we find that the downconversion spectrum 
is limited by the dispersion across the signal freqüèn- 
cies. This differs from the example studied in Ref. |22T | 
which assumed noninstantaneous second-order nonlinear- 
ity leading to limited downconversion spectrum, while 
the dispersion within the bandwidth of down-converted 
freqüències was neglected. Also, we discuss in detail the 
perturbative regime of weak conversion efficiency which 
provides intuitive insights into the structure of squeezing 
building on previous results obtained for spectral two- 
photon wave functions [2M E3] . We also show that 
the characteristic squeezing modes of a parametric ampli- 
fier provide a general answer to the question of optimiz- 
ing the local oscillator shape in homodyne measurements 

HI EU. 

This paper is organized as follows. In Sec. |n] we give 
theoretical foundations of the decomposition used to ana- 
lyze the structure of squeezing in the parametric process. 
Then in Sec. lIIII we discuss the weak pumping regime with 
a low pair creation probability per pulse and show that 



the squeezing modes are defined primarily by the mate- 
rial dispersion coefficients. The general case of strong 
pump fields which needs to be solved numerically is ad- 
dressed in Sec. IIVI The consequences of the multimode 
character of the output field for homodyne detection are 
studied in detail in Sec.0 Finally, Sec. lVIl concludes the 
paper. 



II. GENERAL MODEL OF AN OPA 

The subject of our analysis is one-dimensional propa- 
gation of an optical field through a nonlinear medium in 
which three-wave mixing takes place. For simplicity and 
specificity, we will focus here on the one-dimensional case 
that can be realized in non-linear waveguides [sï}, which 
recently have become available commercially. The one- 
dimensional case nevertheless includes all the complexi- 
ties of multimode propagation in the spectral degree of 
freedom, and the generalization of the formalism intro- 
duced here to the full three-dimensional case is straight- 
forward. 

If the pump field is strong enough to neglect its deple- 
tion as well as fluctuations, the propagation equations in 
the Heisenberg picture become linear in the signal field. 
Following the Standard nonlinear òptics description of the 
evolution of the fields with respect to the distance covered 
rather than time, we shall adopt the modal decomposi- 
tion of the elèctric field operator in the form 



, duj 

E(z,t)=i / — 



TlLÜ 



2ttV 2e ^R 



e- luJt à{uj; z) + H.c. (2) 



where eo is the permittivity of the vacuum, n(u>) is the 
refractive index of the medium at a given frequency and 
polarization direction, à(u; z) is the annihilation oper- 
ator of a monochromatic mode at a given propagation 
stage indicated by z. We have neglected the dependence 
of the field on transverse coordinates, as we concentrate 
here on the one-dimensional case of a waveguide. The 
quantization of the classical propagation equations of the 
three-wave mixing problem 3.1] is done by a formal sub- 
stitution of classical amplitudes by field operators. This 
yields a first-order differential equation for annihilation 
operators à(u>; z) of the form: 



dà(uj; z) 
dz 



= ik(uj)à(u>; z) 



(3) 



+— i— / du 1 é k * {u >' + ^ z E p {u' + Lü)à\u';z) 

where the first term on the right-hand side represents 
the free propagation of the field in the medium, while the 
second one represents nonlinear interaction. In the above 
expression, E p (u>) is the spectral amplitude of the pump 
field at the point z = 0, which is chosen half-way through 
the waveguide of a total length L, ui p is the central fre- 
quency of the pump spectrum, k(u>) and k p (u>) are the 
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signal and the pump field wavevectors, respectively, and 
we rescaled the spectral amplitude of the pump field by 
introducing E — J diüE p (üj). For all calculations in this 
paper, E p (w) is assumed real and positive, i.e. the pump 
pulse is prepared so that in the middle of the waveguide 
it is not chirped. 

With this choice, the multiplicative constant Lnl is the 
characteristic length of the nonlinear interaction defined 
as 



1 



L NL 8c 2 fc(ay 2) ' 



(4) 



where d e g is the effective nonlinearity coefficient for the 
given polarization directions in the medium and io p /2 is 
assumed to be the central frequency of the signal field. 
When introducing Lnl, we assumed that the frequency 
dependence of d e g combined with the normalization fac- 
tors in the modal decomposition given in Eq. yields 
a constant whose changes over the bandwidth of the sig- 
nal field can be neglected. The physical interpretation of 
Lnl is that when the signal field is restricted to a single 
frequency mode and no phase mismatch with the pump 
field is present, Lnl is the length over which the mini- 
mum and the maximum RMS quadrature noise scales by 
e. This corresponds to the single-mode degree of squeez- 
ing introduced in Eq. (QJ equal to Ç — L/Lnl- Let us 
note that 1/Ínl is proportional to the pump pulse am- 
plitude, which we will further refer to as the pumping 
strength. 

It is instructive to recast Eq. (01 into a Standard form 
which allows us to compare different quantum mechani- 
cal pictures of the evolution of a system. Explicitly, the 
development of the field operators described by Eq. <EJ 
can be written as: 



dà(u); z) 
dz 



= í[H, à(uj; z)] 



(5) 



where the operator H, formally equivalent to the Hamil- 
tonian governing the evolution of the system, is given 
by a sum of two terms H = H + V(z). The first term 
generates linear propagation: 



H = 



du> fc(w)a^ (w)à(w) 



(6) 



while the second term is responsible for the three-wave 
mixing process: 



V{z) 



1 



dujduj' 



2LnlEo 

x [ie tk ^' +uj)z E p (cü' + wjfitfwjfit^') + H.< 



(7) 



The evolution of the quantum state of light while prop- 
agating through the crystal can in principle be evalu- 
ated in the interaction picture by integrating over z the 
term V(z). This is however highly complicated because 
the Hamiltonian is position-dependent and generally 
does not commute with itself at different vàlues of z. 



Consequently, the evolution operator needs to be appro- 
priately ordered in its expansion in V(z). As we discuss 
in the next section and in the Appendix, this ordering can 
sometimes be neglected and decomposition into squeez- 
ing eigenmodes can be performed efficiently using the 
interaction picture. The general case is however best de- 
scribed using the Heisenberg picture approach. 

As the propagation equation given in Eq. <EJ is linear 
in the field operators, for any given waveguide length L 
we can find Green functions C(ui,ui') and S(üj,üj') that 
transform the input field operators a ln (uj) = a(uj; —L/2) 
into output field operators a° ut (uj) = a(uj; L/2) according 
to: 

à out {u) = í dLü'[C(Lü,Lü , )à m (üj') + s(üj,üj')à m \üj')}. 



(8) 

The above Bogoliubov transformation can be brought 
into a canonical form by applying the Bloch-Messiah re- 
duction pïj which we will now briefly review. The rele- 
vant mathematical tool is the singular value decomposi- 
tion applied to the Green functions C(uj, to') and S(lu, ui'). 
Because the output annihilation operators must satisfy 
canonical commutation relations, the functions C(lu,lü') 
and S(lü,uj') are connected through a number of rela- 
tions. In particular, the singular value decompositions 
of both C(lu,lu') and S(lu,lu') can be represented using 
a shared set of parameters and functions. The explicit 
formulas have the following form: 

oo 

C(üJ,U)') = ^ l{>n{u>) CQShÇn (j) n (üj') 

n=0 

oo 

S{u>,J) - ^CHsinhC„CK) (9) 

where 4> n (uj) and ip n (u) are two orthonormal bases and £„ 
are real nonnegative parameters. The existence of such 
a joint decomposition allows us to introduce two sets of 
operators for input and output modes, defined according 
to 

C = f düüM^)à m (tü) 
//,','" = / du,Mu)à out (Lü), (10) 



where ò" 1 and &° uí are the corresponding annihilation 
operators which satisfy the Standard commutation rela- 
tions. Consequently, we can recast the general transfor- 
mation given in Eq. (jHJ into a simple form of single-mode 
squeezing transformations acting in parallel on orthogo- 
nal modes: 



COShCnC+SinhCnC 1 , 



(11) 



bringing us back to the elementary expression given in 
Eq. GJ. 

Since any linear transformation of the form given in 
Eq. can be decomposed into orthogonal modes which 
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OPA 



FIG. 1: Pictorial representation of the Bloch-Messiah theo- 
rem. Any Bogoliubov transformation of input modes à l " into 
output modes à™ can be decomposed into a passive trans- 
formation U — 4>n from modes ó™ into b^™ followed by inde- 
pendent squeezers with gains G„ and another passive trans- 
formation V = ipn from b°n l into à^ ut . 



undergo independent evolution, the above procedure, 
summarized graphically in Fig. can be applied to any 
combination of processes such as linear òptics, optical 
ampliflcation, and Kerr squeezing. Were we able to sep- 
arate modes i/} n (u)) spatially, we could turn a realistic 
parametric amplifier into a generator of múltiple pure 
squeezed states 

Let us note that the procedure of singular value decom- 
position applied above is formally equivalent to Schmidt 
decomposition of a quantum state of a bipartite system. 
In the next section and in the Appendix, this equiva- 
lence will allow us to connect the general decomposition 
to previous treatments of parametric down-conversion in 
the perturbative regime. 

One interesting feature of the squeezing eigenmodes is 
that for real E p (i.e. the pump pulse in the middle of the 
waveguide is not chirped or its chirp is insignificant), 



ip n (uj) = <j)* n {u). 



(12) 



This is proven by noticing the symmetry of Eq. <EJ with 
respect to the point z = 0. Indeed, it can be rewritten 
as follows: 



dà(u>; z) 



-ik(Lü)à(Lü; z) 



(13) 



ing to Eq. 

oo 

C(u,u>') = ^2 ^n{uj) coshC„0*( w ') 

n=0 

oo 

S(u,u') = ^Vn(w) sinh(-C„)^„(w'), (15) 



n=0 



we find, in addition to 110L another set of propagation 
eigenmodes 



yüut 



dLÜ^* n (cü)à m (Lü) 



(16) 



with characteristic squeezing parameters ( n . Assuming 
that all Cn are unequal (which is typically the case as we 
show below), the Bloch-Messiah decomposition is unique, 

and we conclude that b™ — b™ and 6°"' = ò° uí , which 
leads to Eq. 1121 . In the time domain, this identity cor- 
responds to the output modes being the time reversal of 
the input. 

Finally, let us point out that the decomposition given 
in Eq. © has important implications when considering a 
parametric amplifier seeded with a coherent pulse. The 
reason for seeding the parametric amplifier is to impose 
a mode in which the majority of the output photons is 
generated, thus overwhelming emission to other modes 
whose contents can be neglected as an unimportant back- 
ground j^. Eq. lfTT|l indicates that the output of the 
parametric amplifier will contain a pure squeezed coher- 
ent state in a single mode if the seeding coherent state oc- 
cupies one of the characteristic input modes </>„(w). If the 
seeding pulse is in superpositon of many modes <fi n (uj), 
under the action of the parametric amplifier it will gener- 
ate múltiple squeezed coherent modes, in principle each 
one of different intensity and squeezing parameter. The 
decomposition in Eq. Q provides one with characteris- 
tic pulse profiles which are optimal for seeding the 
amplifier. 



which can be formally interpreted as reverse propagation 
of the quantum field from —z = —L/2 to — z = L/2. 
Integrating Eq. l(T3"jl between these two points we write, 
by analogy with Eq. 0: 

à m (uj) = J du)'[C(üj,üj')à° ut {üj) + S(uj,uj')^ out (ij)}, 

(14) 

where the superscripts in and out refer, as previously 
defined, to the points z = —L/2 and z = L/2, re- 
spectively. Comparing Eqs. Q and ll.'H . we establish 
a simple relation between the Green functions of di- 
rect and reverse propagation: C(lü,lü') = C*(üj,uj') and 
S(lü,lü') = —S*(lü,lü'). Decomposing the latter accord- 



III. SINGLE PAIR GENERATION REGIME 

In a general case, the equations of motion for the field 
operators specified by Eq. can be solved only by nu- 
merical means. Before we resort to numerical methods, 
we will discuss in this section an approximate solution to 
the problem in the weak conversion limit, when only the 
first order of the perturbation theory is relevant. This is 
indeed the case if the nonlinear interaction length Lnl is 
much greater that the crystal length L, thus keeping the 
squeezing weak. This approach will give us intuitive in- 
sights into the structure of the output field, and will help 
us to identify parameters relevant to its characterization 
in the arbitrary case. 
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The first order approximation to Eq. © can be ob- 
tained most easily by substituting 



à(u>; z) — exp(ik(uj)z) ài(ui] z) 



(17) 



into Eq. J2l, then integrating it over z and retaining 
terms up to first order in L/Lnl- This procedure yields 
an expression for à{ui;L). Using this result we find the 
approximate Green functions in the form: 



C(lü,lü') = S{uj - Lü')e lk{u)L 



where 



Afc = k p (uj + lü') - k(oj) - k(ui') 



(18) 
LAk 
' 2 

(19) 



(20) 



is the phase mismatch between the pump and the pair of 
down-converted photons. 

In most cases it is sufficient to expand Ak up to the 
second order in deviations from the respective central 
freqüències. Then the approximate expression for the 
phase mismatch takes the form: 



Afc 



(fli# ~ (3i){lü + lü' - ui p ) 

-ift ((w - <ü p /2) 2 + (J - lü p /2) 2 ) 

+ ^2,p{u> + U)' - ^p) 2 , 



where 



0n = 



d n k{uj) 



diü n 

d n k p (Lü) 



LÜ—Lü p /2 



dui r 



(21) 

(22a) 
(22b) 



are dispersion coefficients for the pump and the signal 
fields. The coefncients relevant to our analysis are: the 
difference of inverse group velocities /3i íP — /3i and the 
group velocity dispersion coefficients (3% and p2, P - The 
phase term e l [H^)-H^')\L/2 - m arises form the 

free propagation of the squeezed field in the waveguide. 

Physically, performing the expansion up to the first 
order in L/Xnl means considering only generation of 
single photon pairs. The spectral properties of single 
photon pairs generated by ultrashort pulses have been 
extensively studied in the context of quantum indis- 
tinguishabilitv in two-photon interference experiments 
EaSmEl· The bàsic entity in those studies is the 
two-photon wave function *íf(ui,ui') = l L y|í' ont ) de- 
scribing the probability amplitude of generating a pair of 
photons with freqüències u> and u>'. In the Appendix, we 
relate this wave function to the Green function S(lü,lü') 
through a simple formula: 



0) 




FIG. 2: Gaussian approximation of the Green function 
S(lü,lü'), given in Eq. HlHfl by a product of the pump spectral 
amplitude and the phasematching function. The first factor is 
nonzero on a wide strip centered around the lüs+ojí = lü p line. 
The second factor is nonzero on a strip lying along the Ak = 
line, which is approximately parabòlic and tangent to the line 
lüs + lüí — Lü p at the point (lü p /2,u p /2). Thus the product 
S(iü,iu') is nonzero on a strip ariound the u) s + lüí = u p line 
and extending to points where the phasematching function 
goes off the pump spectral amplitude. 



In other words, the eigenmodes of pulsed squeezing are 
closely related to the eigenmodes of the Schmidt decom- 
position of a biphoton spectrum originally found by Law 
et al. |27j . 

Identity i'2?>\) suggests a Gaussian approximation simi- 
lar to that applied in the analysis of the two-photon wave 
function , which will yield a singular value decomposi- 
tion in a closed analytical form. For this purpose, let us 
assume that the quadratic form of the Gaussian which 
approximates Eq. ifïflll has the principal axes given by 



UJ — Lü 

the form: 



and lü + lü' — Lüp, and seek an expression of 



S< 



V 7TÒA 



x exp 



(lü + Lü' — Lü p ) 2 (Lü — Lü') 2 



2S 2 



2A 2 



(24) 



In order to find the unknown coefficients in the above 
formula, let us take a Gaussian pump pulse of duration 

Tp. 



E p (u>) oc exp 



■^-{lü- Lüp) 2 



(25) 



We insert the expansion given in Eq. (12H into Eq. (I19II . 
apply the following approximations to the sine function 
along the axes ui — lü' = and lü + lü' = lü p : 



9{ui,ui') = S(ui,ui')e ik ^ L . 



(23) 



~ x2 l\ sinc(x 2 ) 



2 /3 



(26) 
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and compare Eq. i|19H after all these approximations with 
Eq. j2Ü) along the principal axes. This yields: 



1 




p ~ 




1 


1 


A2 " 


12 


AT = 


/ 



T 2 



A,p) ! 



Sq{üJ, Ul')\ 2 dLüdLd' 



I? 



4L 2 p 

NL 



2 ÓA 



(27a) 
(27b) 
(27c) 



The meaning of these parameters is depicted in Fig. 
Typically |20], 8 -C A and consequently the bandwidth 
of the downconverted light is approximately proportional 
to A. In the particular case of type I phasematching 
this bandwidth is limited mostly by the group velocity 
dispersion (32 at signal freqüències. Note that, in con- 
trast, the difference between the group velocities at the 
pump and the signal freqüències (3\ — (3\ íP sets no limit on 
the downconversion spectrum. The second relevant fre- 
quency scale is set by the parameter 8 which characterizes 
the width of the spectral correlation between photons in 
pair. This width is limited by the pump pulse bandwidth 
and the group velocity difference (3\ — (3\, v . Also, let us 
note that the Green function S(uj,uj') is not normalized, 
and its contribution is implicitly assumed to be small, of 
the order L/Lnl, as compared to that of C(üj,üj'). The 
normalization constant N of Sg(u, uj') given in Eq. Ü27c|) 
is equal to the total number of photons per pulse emitted 
by the downconversion source. 

Since in general A and 8 are different, one cannot sim- 
ply factorize S(ui, ui') and define a single mode into which 
photon pairs are emitted. If the signal and idler pho- 
tons were nondegenerate, one would call such a spectrum 
entangled; this definition is however inapplicable to the 
present situation because the two fields cannot be sepa- 
rated using another degree of freedom. With the Gaus- 
sian approximation in hand, let us now employ the sum- 
mation formula used previously to find the Schmidt de- 
composition of a bipartite Gaussian state |Ü3,ÜÏ3' which 
can be applied directly to the Green function l|2"3j) in order 
to find its singular value decomposition. Recasting the 
decomposition according to the general expression spec- 
ified in Eq. ^ , we obtain that the squeezing parameters 
Cn form a geomètric sequence given by: 



Ç„ « sinh ( n 



N 



tanh" r 



(28) 



coshr 

with the parameter r = ln(A/<5)/2, while the correspond 
ing eigenmodes are described by the Hermite functions: 



-ik(ui)L/2 



e ik(ui)L/2 
H n (L)T s ) 

\Jl n n\ 



<f>n(u) 



(29) 



-t s 2 (w-w p /2) 2 /2 



with characteristic spectral width of all squeezing modes 
being proportional to the inverse geomètric mean of the 
two spectral widths 8 and A: 



(30) 



These expressions for the eigenmodes and the squeezing 
parameters will serve as a reference when analyzing nu- 
merical results in the nonperturbative regime presented 
in the next section. 

As expected, ipn(u) — 0* (oj), and the modes' spectral 
phases can be attributed to the linear dispersive propaga- 
tion through the BBO crystal. The nonzero phase profiles 
are difficult to handle experimentally. One might think 
they can be made negligible by choosing a sufficiently 
long pump pulse so its dispersion can be neglected. This 
is however not the case because the eigenmodes' charac- 
teristic spectral width 1/t s is determined not only by 8, 
but also by A. The magnitude A is, in turn, determined 
by the properties of the crystal rather than the pump 
pulse and is usually large. 

The fact that the phase profiles of squeezing eigen- 
modes depend on the crystal length has an interesting 
consequence. Suppose we prepare a particular state of 
light in a squeezing eigenmode b™(L) of a crystal of 
length L and insert it into the crystal. After propagation 
through the entire crystal, this state will undergo squeez- 
ing and leave the crystal in the mode ò°"*(L). Now we 
ask, what will be the state of light at the point z = 
half way through the crystal? Propagation of light to 
this point is equivalent to that through a crystal of length 
L/2. The mode b™(L) is however not an eigenmode as- 
sociated with this crystal; therefore, the state of light at 
z = cannot be defined in any pure mode, but only as 
an entangled state of many modes. We see that while 
light enters and leaves the crystal in a pure optical mode 
and state, inside the crystal all input modes temporarily 
become entangled with each other. 



IV. INTENSE GENERATION REGIME 

In this section, we discuss numerical solutions to the 
propagation equation in a general, not necessarily 
perturbative regime. As a concrete example, we will 
consider type-I interaction in a BBO crystal. We as- 
sume that a waveguide is formed in a direction of perfect 
phase matching for an interaction of a 400 nm pump wave 
with 800 nm signal waves (corresponding to the angle be- 
tween the waveguide axis and the optical axis 9 = 29.2°). 
For this system, we have solved numerically Eq. © with 
k(uj) and k p {u>) approximated by vàlues for a bulk BBO 
crystal, and a Gaussian pump field given by Eq. 
Let us stress that in the calculations presented below we 
use neither the crude Gaussian approximation of the sine 
phasematching function, nor are we limited to the first 
order aceuracy in L/Ínl- Nevertheless, the simple model 
developed in the preceding section will guide us through 
general solutions. 

As the propagation equation is linear in the strong 
undepleted pump approximation, Eq. <J2J is identical to 
the classical equation of motion for a pulse propagat- 
ing through a waveguide, with the annihilation opera- 
tors à(uj; z) replaced by the spectral amplitudes of the 
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elèctric field at respective signal freqüències a(uj;z) [3Ï1 |. 
Thus we can compute numerically the matrix approxima- 
tions to the Green functions C(ui,ui') and S(oú,üú') using 
the well-established in classical nonlinear òptics split-step 
method for solving partial differential equations. The 
approximations are obtained by solving the propagation 
equation for a complete set of <5-function shaped initial 
conditions where the classical signal field a(u>; z = —L/2) 
is taken to be equal to zero everywhere on the computa- 
tional grid in the frequency domain, except for a single 
point of the grid at which it assumes 1 or i. When the 
distinguished point corresponds to the frequency uj 1 , solv- 
ing the propagation equation with such initial conditions 
yields single columns of the discretized Green functions 
C(ui,ui') and S(ui, ui'). After calculating all the columns 
of the Green functions, we computed the singular-value 
decompositions defined in Eq. © which gives the desired 
Bloch-Messiah reduction. 

We assessed whether the computational grid was fine 
enough to support all the modes with nonunit gain by 
checking that the squeezing parameters for high-order 
modes approach zero. Additional tests of computational 
accuracy consisted in verifying that the singular vàlues 
of C(w, u/) and S(u>, u>') are pairwise sinh and cosh func- 
tions of the same real number, as stated in Eq. ^ and 
that the calculated modes uphold the relation ltÏ2ll . Al- 
though many modes are not included in our numerical 
calculation, we argue that their squeezing is negligible, 
thus the actual choice of the mode functions (j> n (to) and 
ip n (ijj) representing them is unimportant. 

Following the approach used in nonlinear òptics, the 
calculations were performed in the reference frame of 
the down-converted light, moving with the group velocity 
l/Pi- This corresponds to substituting 

Ò(w; z) -» e M°+Pi("-"p/V]* à{uj; z) (31) 

into the propagation equation given in Eq. <EJ| . The mode 
functions <j> n (u)) and ifj n {uj) will be given also in the mov- 
ing reference frame. This reduces to stripping the mode 
functions from the linear phase in the expansion around 
ujp/2. Experimentally, the linear phase can be always 
corrected by adjusting the temporal delay of the down- 
conversion beam, whereas higher order phase terms lead 
to pulse distortions requiring dispersion management. 

In our numerical calculations we used a Gaussian pump 
pulse of length t p = 24 fs (40 fs intensity FWHM) and the 
waveguide length L — 1 mm. This corresponds to the in- 
verse bandwidth of the downconversion light A" 1 ~ 4 f s 
and the inverse of the spectral width of correlations be- 
tween the two photon in a pair <5 _1 ~ 50 fs. We will study 
how the features of the down-converted light depend on 
the nonlinear interaction length Lnl- 

We found that up to L/Lnl = 15 (which corre- 
sponds to over 100 dB maximal quadrature squeezing) 
the squeezing parameter of any squeezed mode is in- 
versely proportional to nonlinear length with a specific 
proportionality constant. These vàlues are depicted in 
Fig. 03 We note a substantial discrepancy with the Gaus- 
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n 

FIG. 3: The rescaled squeezing parameters LnlCïi as a func- 
tion of the mode number, verified numerically to be indepen- 
dent of Lnl for L/Z/NL < 15. Empty bars depict results of 
the numerical computation, while black bars represent pre- 
dictions of the Gaussian approximation. 



sian model from Sec. lIIIl which predicts that the squeez- 
ing parameters should follow a geomètric sequence (l28l) . 
represented by a straight line in the semi-log scale plot 
in Fig. □ 

The Gaussian approximation turns out to be much 
more successful in predicting the shapes of the charac- 
teristic squeezing modes. Particular ly in the weak pump- 
ing regime, the numerically calculated input and output 
modes (f> n (u>) and ip n (w) f° r l° w n are almost identical to 
the Hermite-Gaussian functions predicted by Eq. (l29l) . 
with only a slight asymmetry and more abrupt decay 
of their outermost wings (Fig. QJ • For strong pumping, 
when the nonlinear interaction length Lnl becomes com- 
parable or smaller than the crystal length L, the spectral 
intensity profiles begin to broaden and change in shape 
(Fig.[S]| due to nonlinear modulation of the optical fields. 
Still, within the range of the interaction lengths studied 
here (i.e. L/Lnl < 15), the Hermite-Gaussian approx- 
imation for the several most strongly squeezed modes 
works reasonably well. We expect this to be the case 
as long as the higher order terms in the phase mismatch 
expansion l(2Ï|l are negligible. In the case of BBO, this 
condition is satisfied for t p ,t s 3> 10 fs. 

In order to characterize the change in the mode 
widths, we fitted the numerically obtained spectral in- 
tensities (t^)| 2 with the Gaussian-Hermite approxima- 
tions Eq. i|29ll for each Lnl and n, treating the character- 
istic time t s as a free parameter. We found this param- 
eter to deviate from the constant 20 fs value predicted 
by Eq. ll-'iOL and, as evidenced by Fig. El to depend on 
both the mode number and the pumping strength. In the 
weak pumping regime, t s « 15 fs for n = and reduces 
by about 3% for each subsequent mode, and is indepen- 
dent from 1/nl- For stronger pumping (Lnl fí 1), t s falls 
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FIG. 4: Spectral intensities |7/) n (o;)| 2 of the dominant three 
modes n = (sòlid line), n = 1 (dashed line), and n = 2 
(dotted line) for Lnl = 100 mm. 




frequency [THz] 

FIG. 5: Spectral intensity and phase of the maximally 
squeezed mode i/jq(uj) for Lnl = 100 mm (sòlid and dashed- 
dotted line respectively) and Lnl = 1/15 mm (dashed and 
dotted lines). For the case studied, the input modes are com- 
plex conjugate to output modes, </>o(lü) = $J(w). 



quickly with the increasing pumping strength. 

As expected, we found the spectral phases of the in- 
put and the output modes to have opposite signs with 
virtually equal absolute vàlues. For strong pumping, the 
phase profile also becomes flatter for shorter interaction 
lengths Lnl, as illustrated in Fig.0 Thus in the intense 
generation regime we can no longer attribute the spectral 
phase of the modes to linear dispersive propagation. All 
these effects mean that the features of multimode squeez- 
ing become more fragile compared to the single-pair gen- 
eration case. Thus the exact modal characteristics need 
to be carefully taken into account when manipulating or 
detecting strongly squeezed light. 



16 : 




FIG. 6: Duration of the dominant squeezed modes r Sj o (sòlid 
line), t s ,i (dashed line) and r Sj 2 (dotted line) as a function 
of I/Lnl, obtained by fitting Hermite-Gauss profiles from 
Eq. I|29jl to numerical results. 

V. HOMODYNE DETECTION 

Homodyne detection is a fundamental diagnòstic 
and measurement technique in applications of squeezed 
states for quantum-enhanced metrology and continuous- 
variable quantum information processing. It is a difficult 
measurement because of the electrònic noise, losses in the 
paths of the squeezed beams and nonideal mode matching 
between squeezed field and local oscillator which defines 
the measured quadrature. This section is devoted to a 
detailed study of the latter issue, and the modal decom- 
position of a parametric amplifier introduced in Sec. |n] 
provides a natural framework for its discussion. 

We will consider an unseeded optical parametric am- 
plifier followed by a balanced homodyne detector with 
the local oscillator (LO) pulse prepared in a certain 
shape characterized by the mode function i/ilo( w )- The 
quantity of interest is the measured quadrature noise. 
The analysis is straightforward in the particular case 
when the local oscillator employed in the homodyne mea- 
surement is prepared in one of the characteristic modes 
V'lo(w) = ip n (uj). This choice of the LO means that we 
are "listening" only to one of the independently squeezed 
characteristic modes of the amplifier, whose evolution is 
given by Eq. itTTll . For vacuum input, this mode evolves 
into a squeezed vacuum state with the squeezing param- 
eter equal to £ n . Thus the detected mean square quadra- 
ture noise (Q 2 s ) as a function of the LO phase 9 is equal 
to Q: 

(Q 2 n 9 ) = I( e 2 í" sin 2 9 + e- 2Ç " cos 2 9). (32) 

In particular, the product of the maximum (Q\) = 
e 2Ç "/4 and the minimum (Q'L) = e _2 ^™/4 quadra- 
ture noise will be the same as for the vacuum state 
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{Q\){Q 2 -) = 1/16 — the minimum value allowed by the 
Heisenberg uncertainty relation. 

In the general case of an arbitrary LO pulse shape we 
can calculate the quadrature noise level by decomposing 
V'lo(w) in the basis of the characteristic output modes 
of the amplifier: 

oo 

Í>LO (W) = H , (33) 

n=0 

where M n and 9 n are real numbers and we have as- 
sumed that the spectral amplitude ipLo{u) is normalized 
to unity. It is important that the quantum quadrature 
fluctuations of different characteristic modes are inde- 
pendent because the mutually uncorrelated modes b™ 
independently evolve into 6° uí . Thus, using the local os- 
cillator of the form l-i-'tll , we measure a sum of noise con- 
tributions from the quadratures (Q^ g ) weighted with 
Aí 2 - 

00 M 2 

(Qlo) = E -f ( e2C " sin2 e ™ + e ~ 2C " cos2 «) ( 34 ) 

n=0 

In general, the above expression combines contributions 
from both squeezed and antisqueezed quadratures. How- 
ever, if all the phases 9 n are equal, we can always bring 
them to by introducing an additional phase delay in 
the LO arm. Then the minimum quadrature noise that 
can be observed is below the shot noise limit: 

00 M 2 1 
(Q-) = E^ ex P(- 2 ^)^4- ( 35 ) 

n=0 

The condition of equal phases 9 n can be easily realized in 
the single pair generation regime, because all the output 
modes Eq. (129)) bear the same spectral phase associated 
with the linear propagation in the waveguide. Thus for 
an unchirped local oscillator with an arbitrary intensity 
spectrum Ilo{u) we can optimize the amount of detected 
squeezing by making its phase equal to the phase ac- 
quired by the squeezed field during propagation through 
a crystal of length L/2: 

^lo (w) = y/holü) e^ L ' 2 . (36) 

Let us note that if some of the squeezing parameters 
C, n are equal — for concreteness let us take £o and Ç,i 
— then any local oscillator which is in a phase locked 
superposition of the corresponding modes tpLo(v) — 
MQipo(u})+Miipx(uj) with real Mq and M\ will enable the 
detection of the same minimum quadrature noise equal 
to (Q 2 _) = exp(-2Co)/4 = exp(-2Çi)/4. For the weak 
pumping regime, several largest Ç n are indeed approxi- 
mately (albeit not exactly) equal, as shown in Fig. 03 As 
a result, we have a substantial freedom in choosing the 
the local oscillator mode among various linear combina- 
tions of primary squeezing modes without affecting the 
detection efHciency. 



When the squeezing becomes more intense, the choice 
of the local oscillator mode becomes more and more criti- 
cal in order to detect squeezing. This is easily seen in the 
Standard single- mode case described by Eq. i'.VHi . as even 
a small deviation from 9 — adds a significant amount 
of noise from the antisqueezed quadrature. In the case 
of the local oscillator distributed over several characteris- 
tic output modes of the amplifier, detection of squeezing 
requires accurate control of each of the terms in its de- 
compositon (l33f) . This is illustrated with Fig.Q depicting 
the observed squeezing for various nonlinear interaction 
lengths Lnl, based on the numerical results of section IIVI 
The squeezing is calculated as a function of the inverse 
bandwidth tlo of the local oscillator pulse, assumed to 
have a Gaussian spectral intensity profile with a phase 
satisfying Eq. I|36JI : 

toM - ^exp (-^o+iL·k^ . (37 ) 

It is seen that in the single-pair generation regime, when 
i/ííNL <C 1, the squeezing is detected for a broad range 
of local oscillator lengths. The local oscillator duration 
tlo can thus be swept approximately from the inverse 
bandwidth of the downconversion light A _1 ~ 4 fs to 
the inverse of the spectral width of correlations between 
the two photon in a pair <5 -1 ~ 50 fs. 

This result permits a physical interpretation of down- 
conversion in the time-domain picture. The two photons 
in a pair can be born at any moment in time within the 
duration of the pump pulse "smeared" by a difference 
of the group velocities between pump and signal pulses 
Pi.p — fil [see Eq. p7a| ]„ i.e. within the time interval 
1/6. If the LO pulse of a longer duration is chosen, it 
is known a priori to contain modes into which no pho- 
tons have been emitted, thus entailing inefficient detec- 
tion. The existence of the upper limit A to the local 
oscillator bandwidth is explained by the finite length of 
the down-conversion crystal. Due to the group velocity 
dispersion (which limits the down-conversion spectrum), 
single-photon wavepackets propagating through the crys- 
tal di verge in time by A _1 . If the local oscillator pulse 
is chosen too short, it may happen that one photon in 
a pair is registered within the local oscillator mode, but 
the other one arrives either sooner or later, also leading 
to a reduction in detection efficiency. 

As seen in Fig. when the pumping intensity in- 
creases, reducing Lnl, the choice of the LO duration tlo 
becomes more and more critical and finally, due to the 
mismatched spectral phase, no choice of tlo matches the 
most strongly squeezed characteristic mode. Moreover, 
it turns out that using an LO pulse significantly longer 
than the fundamental squeezed mode is favorable. This 
effect appears paradoxical, because the squeezing modes' 
characteristic time in fact reduces with stronger pumping 
(Fig. El, but it can be explained as follows. Numerical 
calculations show that for t L o — T s the decomposition 
il33)l contains about a 2% contribution from modes with 
odd n, as shown in Fig. |S| We found that for small n, 
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FIG. 7: Minimum quadrature noise detected using the lo- 
cal oscillator given by Eq. 13711 of duration tlo for nonlinear 
lengths Lnl = 1 mm (dotted line), Lnl = 1/2 mm (dashed- 
dotted line), Lnl = 1/3 mm (sòlid line), and Lnl = 1/4 mm 
(dashed line). 




FIG. 8: Squared coemcients M% for the decomposition of 
a local oscillator pulse given by Eq. I|37fl into characteris- 
tic squeezing modes according to Eq. \A',ty . for the nonlin- 
ear length Lnl = 1/2 mm and tlo = 15 fs (black bars), 
tlo = 30 fs (grey bars), and tlo = 50 fs (empty bars). Note 
the logarithmic scale on the vertical axis. 



corresponding to the most strongly squeezed modes, the 
phases in the decomposition in Eq. follow approx- 
imately the rule 6 n ~ mr/2. Consequently the modes 
with odd n, despite having only a minor contribution to 
LO shape, build up altogether significant noise, as we are 
observing their noisy quadratures, described in Eq. l(34*jl 
by terms proportional to . For tlo — 2t s , on the other 
hand, even though many more modes contribute to the 
LO decomposition in Eq. 133[l . they are predominantly 
even with equal phases 8 n ~ 0. According to our simu- 
lations, the loss of squeezing at high pumping intensities 
is present vàlid even if we allow more general LO pulse 
shapes with a Gaussian spectral profile and an arbitrary 
quadratic phase. 

The quadrature noise observed in homodyne detection 
of multimode squeezed light is similar to that of a single- 
mode squeezed state with an inefficient detector in that 
it does not exhibit minimum uncertainty: {Q\){Q 2 _) > 
1/16. In the latter case, the efficiency of the detector can 
be calculated as follows: 



V = 



-16(Q^)(Q 2 _)+4(Q^)+4(Q 2 _)-1 



0.98 



0.96 




0.94 



0.92 



(38) 



FIG. 9: Quantum efficiency r\ of detecting squeezing defined 
in Eq. Q38fl as a function of the local oscillator duration tlo 
for nonlinear lengths Lnl = 1 mm (sòlid line), Lnl = 1/2 mm 
(dashed line), Lnl = 1/3 mm (dashed-dotted line), and 
Lnl = 1/4 mm (dotted line). The thin dashed line repre- 
sents the result of the Gaussian approximation. 



It is sometimes useful to express the quality of mul- 
timode squeezed light in terms of the efficiency calcu- 
lated according to the above relation. In Fig. El we plot 
this quantity for a local oscillator of the shape given by 
Eq. I|37|) . as a function of the LO pulse duration, for 
several squeezing strengths. We observe the same gen- 
eral behavior as in Fig. In the single-pair generation 
regime, high-77 squeezing detection is possible in a broad 



span of LO modes within the interval 1/A < tlo < 

as the pumping intensity is increased, the choice of LO 

becomes a lot more critical. 

Most existing experiments on pulsed squeezing derive 
the local oscillator from the master làser whose second 
harmònic serves as a pump for the down-conversion. Ac- 
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cording to Figs.0and|HJ this is not always the best strat- 
egy. In the weak pumping regime, the highest efficiency 
is achieved for the LO pulse duration of t s , which is usu- 
ally much shorter than the fundamental pulse width ll-'UH . 
With increasing squeezing strength, the optimal LO pulse 
width increases. 

Combining these considerations allows one to match 
the parameters of the setup, such as the duration of the 
fundamental and second-harmonic pulses, as well as the 
characteristics of the down-conversion crystal, in order to 
reach the optimum detection of squeezing. 



VI. CONCLUSIONS 

In summary, we have analyzed the spectral properties 
of pulsed squeezed light generated by means of single-pass 
optical parametric amplifier. The main tool in our anal- 
ysis was the Bloch-Messiah decomposition of the Green 
functions describing the parametric down-conversion pro- 
cess. This decomposition allowed us to identify the char- 
acteristic squeezing modes, which provide a simple way 
of describing quantum statistical properties of the multi- 
mode output field. 

We have developed an analytical model describing the 
perturbative regime of single pair generation L/Lnl <C 1, 
when the nonlinear interaction length Lnl substantially 
exceeds the crystal length L. The model yielded an ap- 
proximate form of the squeezing modes as well as their 
characteristic time scale. In this regime, high-efnciency 
squeezing expands over a broad range of pulsed modes 
temporally synchronized with the pump pulse. Some of 
the results of the analytical model, such as the mode 
shapes, remain qualitatively vàlid also beyond the per- 
turbative limit. 

For the multiple-pair generation case we have pre- 
sented realistic numerical calculations of the characteris- 
tic modes of a parametric amplifier, in which pure squeez- 
ing is present. Beyond the perturbative regime, the mode 
functions depend on the pumping intensity, including the 
change in their spectral phases. As with increasing para- 
metric gain the fluctuations in antisqueezed quadratures 
grow exponentially, it becomes more difficult to elimi- 
nate their parasitic contribution to the observed homo- 
dyne noise. Therefore, one needs precise tailoring of the 
local oscillator mode, including both its amplitude and 
its phase, if substantial squeezing is to be observed. 

The multimode character of the output field generated 
by an optical parametric amplifier has important con- 
sequences when considering quantum information appli- 
cations based on continuous variables, such as quantum 
cryptography. For example, if a quantum communica- 
tion protocol involves modulating squeezed light, care 
must be taken to ensure that only the mode detected at 
the receiving station experiences modulation. Otherwise 
other modes might carry the same signal and could be 
detected by a third party without causing any observable 
disturbance. Furthermore, if the quadrature fluctuations 



detected by the legitimate user remain above the min- 
imum uncertainty limit, this implies that the observed 
mode was correlated with other undetected modes which 
again could be used to gain information by a third party. 
The multimode nature of pulsed squeezing also needs to 
be taken into account is quantum state engineering com- 
bining discrete- and continuous- variable methods, such 
as preparation of squeezed single-photon states [2fl tóa. 
Schròdinger cats [4£|], and entanglement purification |41| . 
When multimode squeezing is undesired, it can be pre- 
vented by choosing down-conversion configuration with 
unentangled biphoton spectrum |3fiL I42I Uflj . 

Let us also note that the decomposition of the output 
of the parametric process into characteristic squeezing 
modes enables a straightforward analysis of the result 
of combining two squeezed beams on a balanced beam 
splitter. It is well known that in the single-mode picture, 
such a procedure yields a pair of twin beams || . In the 
realistic case, if the two incident squeezed beams are gen- 
erated by identical parametric amplifies, they will com- 
prise the same characteristic squeezing modes perfectly 
matched pair-wise. Consequently, combining them on a 
beam splitter will yield multimode twin beam sets which 
exhibit minimum-uncertainty quadrature correlations in 
the characteristic modes of the squeezers. 
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APPENDIX A: ANALYSIS IN THE 
INTERACTION PICTURE 

The evolution of the quantum state of light while prop- 
agating through the crystal follows: 



|V>r í >=^exp 



L/2 

-i J Vi{z)dz 

-L/2 



W) , (Al) 



where the subscript / denotes the interaction picture and 
T denotes z-ordering, analogous to Standard time or- 
dering. The three-wave mixing Hamiltonian can be 
rewritten in the interaction picture as 



Vj[z) 



1 



duj du 



2Lnl-Eo 

x[íe tAkz E p (uj' + u)à\üj)à f (Lü') +H.. 



(A2) 



where Afc is given by Eq. 1I2Ò")) . In the weak conversion 
limit, when only the first order of the perturbation theory 
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is relevant, we can neglect the z-ordering and obtain 

,^ f dujdu'^i (lj.lj 1 ) [à^ (oj)à^ (üj') — à(Lj)à(uj')] \„LÍn 



\Í>T*) = e 5 
where the integration kernel 



(A3) 



T . LEJlú + u') . LAk 
= p sine— - 

determines the frequeney correlations of photons in a 
down-converted pair. In the weak excitation limit with 
the vacuum input, the output state takes the form: 

|V/ UÍ ) = |0) + \ J dudu/*r(u,u/) \l u , I*/) . (A4) 

In the Schrodinger picture, the output state aequires 
an additional optical phase: 



* s (w,w') = * / (^, W ')e l[ ' cM+fe( "' )li/2 - 



We recover Eq. (GBl from Eqs. JE), l|A4|) . and iXHjl . 

Note that the biphoton spectrum can be diagonalized 
according to Eq. 11411 



* S (w,ü/) = £Cn^n(w)^n(w'). (A6) 

and the output state can be rewritten as follows: 

|^ uí > = |0) + ^C«(C ít ) 2 |0), (A7) 



(A5) 



n=0 



where the b n 's represent individual squeezing eigenmodes 
ltlH)l . As we see, the problem of finding the characreris- 
tic modes of a squeezer in the weak pumping limit is 
equivalent to decomposing the biphoton spectrum into 
uncorrelated components. 
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